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ABSTRACT 



Planetary nebulae (PNe) have diverse morphological shapes, including point- 
symmetric and multipolar structures. Many PNe also have complicated internal 
structures such as torus, lobes, knots, and ansae. A complete accounting of all 
the morphological structures through physical models is difficult. A first step 
toward such an understanding is to derive the true three-dimensional structure 
of the nebulae. In this paper, we show that a multipolar nebula with three 
pairs of lobes can explain many of such features, if orientation and sensitivity 
effects are taken into account. Using only six parameters — the inclination and 
position angles of each pair — we are able to simulate the observed images of 
20 PNe with complex structures. We suggest that the multipolar structure is an 
intrinsic structure of PNe and the statistics of multipolar PNe has been severely 
underestimated in the past. 

Subject headings: planetary nebulae: general 
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Introduction 



One of the fascinating aspects of planetary nebulae (P Ne) is their di verse morphological 
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9181). many attempts 



Stanghellini et al. 



terms such as "round" , "elliptical" , and "bipolar" . As the sensitivity of imaging techniques 
improved, increasing numbers of PNe were found to possess more than one pair of bipolar 



lobes. In early literature, the terms "twofold bipolarity ", "over 



and " bi-butterfly" were used to describe these objects ( jPascoli 



apping bi polar structures" 



1990aU bh 



Manchado et al. 



( 119961 ) introduced a new morphological class named "quadrupolar" : lobe pairs symmetric 



about two distinct axes. In general, the term "multipolar" refers to morph ologies with 



at lea st two axes of symmetry. It is sometimes referred to as "polypolar" ( lLopez et al. 
19981 ). The multipolar morphology is not only found in PNe, but is also observed in earlier 
phases of post-asymptotic giant branch (post-AGB) objects including pre-planetary nebulae 
(pPNe). 

Deeper imaging of PNe has revealed that some well-known PNe are actually multipolar. 



For ex ample, both NGC 2440 and NGC 6072 were classified as bipolar by 



Corradi &: Schwarz 



( 1995 ). However, detailed spatial and kinematical studies of optical lines by 



Lopez et al. 



( 119981 ) confirme d the presen c e of t wo pairs of lobes in NGC 2440, also identified in CO 



observations by 



Wang et al. 



(120081 ). Similarly, NGC 6072 appears to be an elliptical 
shell with faint extensions in the Digital Sky Survey optical image, but the Spitzer Space 
Telescope Infrared Array Camera (IRAC) and Canada— France— Hawaii Telescope (CFHT) 
molecular hydrogen ( H.?) images reveale d a complicated structure with three pairs of lobes 



and an elliptical ring (IKwok et al. 



at least two bipolar outflows ( jMeaburn et al. 



20101) . NGC 729 3 (the Helix Nebula) is known to have 



20081). NGC 6853, the Dumbbell Nebula, 
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shows a pair of cones in addit ion to the well kn own dumbbell-shaped main nebula in the 



IRAC and CFHT H 2 images flKwok et al 



20081 ). In these examples, it can be seen that 



deeper exposures and observations in other wavelengths can result in the identification of 
multipolar lobes. 

Statistical results of the PNe morphology also show an increasing fraction of the 
multipolar class. In early classifications, there were no multipolar categories, and such 
objects might have been included in the "point-symmetric " or "irregular" classes (e.g., 



Balick 



1987 



Schwarz et al. 



1993 



Corradi fc Schwarz 



1995). In more recent results, from 



Hubble Space Telescope (HST) images, 12 % of 140 post-AGB obj ects and 20% of 119 



young PNe were classified as multipolar by 



Manchado et al. 



(120111 ) and 



Sahai et al. 



(2011) 



respectively. The trend may still be growing. 

Efforts have been made to c onstruct three-di mension al (3D) models of indivi dual 



multip olar PNe: e.g. 



J 320 (IHarman et al. 



NGC 6644 flHsia et al. 



20101), Hb 5 (IGarcia-Diaz et al 



201l|), and 



20041 ) . However, as these models are tailored- made for each PN, the 
similarities and differences among objects are unclear, which hinders our understanding of 
the origin of multipolar morphology. In order to explain the presence of multiple outflow 
axes, one has to introduce additional hypotheses, such as precessing jets in the short 
transition period between the spherically symmetric AGB phase and the multipolar post- 



AGB phase. Bipolar, r otating episodic jet s provide dire ct evidence o 



observed in Fleming 1 (ILopez et al 



19931 ) and KjPn 8 (ILopez et al. 



prec ession motions as 



19951 ). It is suggested 



that a rotating jet can account for the creation of continuous point-s ymmetric features at 



different distances fro m the geometric cent er, e.g., IRAS 16585-2145 (j Guerrero et al. 



and IRAS 17028-1004 f Corradi et al. 



maser spots in early stages (e.g. 



20081) 



201 lj) . More over, precessing motions traced by water 



If 



Yung et al.ll201ll ) may bring out an indirect way of proof. 



It is still under debate whether the multiple lobes are formed simultaneously or episodically 
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( ISahail 120021 ). They may involve totally different physical processes. 



Before starting to establish the theories, the first step should be to determine the actual 
3D structure, rather than applying only the projected two-dimensional (2D) images of these 
PNe. Based on the 3D model, one can estimate the kinematic timescale in each outflow 
direction to verify whether they were produced at the same time. It has been pointed out 
that the current classification schemes do not reflect the true 3D physical structure of the 
nebulae as the observed images suffer from effects of sensitivity d ependence, atomic or 



20101 ) . As a first step toward 



molecular species dependence, and projection effects ( iKwokl 
understanding the true 3D structure of PNe, it would be useful to build a unified 3D model 
to see how far one can go to reproduce the observed 2D images of individual objects. By 
making this attempt, we are trying to address the question of whether "each PN is unique" 
or whether they are just different manifestations of a unified structure. 

In this paper, we explore answers to this question by beginning with a simple model of 
a PN with multiple lobes that all have the same shape, kinematic timescales and outflow 
velocities. A sample selection of real PN images to be compared is in the next section, and 
the model is described in detail in Section [3j Section @] presents the projected images of the 
model in different parameters and a comparison with real PN images. Discussions including 
the meanings of parameters and the formation mechanisms can be found in Section [5] 



2. Sample Selection 



To co mpare the modeled images with re a 



literature (iSahai et al. 



2011 



Ueta et al 



2007 



observed ones 


, we 


Guerrero et al. 


2008 



i the 



Harman et al. 



2004 



and references therein) for objects which have been observed in Ha with the Wide-Field 
Planetary Camera 2 of the HST. Emission in the recombination line Ha clearly indicates 
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that the sources are PNe rather than pPNe in which the central stars are not hot 
enough for photoionization. Most of the objects ar e called young PN e according to their 



[O lll]A5007/Ha flux ratio of less than about unity (jSahai et al. 



20111 ) . Ha emission, as the 



most commonly observed line of PNe, provides a larger pool for sample selection in this 
paper and future projects. The sample contains a total of 20 PN images retrieved from the 
Canadian Astronomy Data Center archive (details are listed in Table [2]). These 20 objects 
are selected due to their resemblance to the projected images of our model. Therefore, we 
do not claim that they represent the whole PN set. 



The Model 



Our model images are constructed using the 3D modeling software SHAPE (jSteffen 



20111 ). SHAPE is an interactive tool for morpho- kinematic modeling and image 
reconstruction. Customized objects can be constructed by starting with built-in basic 
geometrical shapes and adding modifiers to them. 

We begin with a simple model where the nebula consists of three pairs of identical 
lobes. We first address the issue of projection effects on the lobes as observed from different 
orientations. All other parameters such as the density profiles are fixed and the inclination 
angle i and position angle (PA) of each pair are the only variable parameters, i is taken to 
be the angle between the lobe axis and the line of sight (LOS). As a result, there are only 
six independent parameters. From these six parameters, the separation angle 9 between 
any two pairs of lobes can be calculated from the inner product. There are three values of 9 
for three pairs of lobes. The lobes are hollow inside with evenly distributed density within 
the "walls" of the lobes. The brightness is proportional to the square of the column density 
based on the fact that the inte nsity of H re combination lines is directly proportional to the 



square of the electron density (IKwok 



20071 ) 
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4. Results 

4.1. Apparent Morphologies 

Figure [1] shows the projected images of 100 randomly generated combinations of 
the three-paired model. Each pair of lobes in each image has an independent random 
orientation in the 3D space: PA is a random multiple of five between —175 and +180 
inclusive (i.e., —175, —170, —165, 0, +180); i is a random integer between and 
90 inclusive with a weight proportional to sinz, as the differential solid angle in polar 
coordinates is (sinz d(i) d(PA)). This expression indicates that a pair of lobes has a higher 
chance to be viewed edge-on than pole-on. The values of i and PA are listed in Table [TJ It 
can be seen that although the true sizes of lobes are the same, the apparent length changes 
with the viewing angle. The combination of three pairs greatly increases the variance. In 
some special cases, the multipolar nature of the nebula is not obvious from the projected 
images, for example, when projections of two or more pairs of lobes are aligned along the 
same direction so that the lobes overlap with each other as seen by the observer. If one 
pair is viewed nearly pole-on or slightly tilted in the equatorial direction, then it may be 
wrongly interpreted as a feature in the equatorial direction or be labeled as a torus. 

From this simulation, we can estimate the fraction of apparent morphologies created 
as the result of random orientation. The visual appearance of these images can be classified 
into the following classes: 

• Class 0: a round or elliptical appearance and no obvious lobes. The fraction belonging 
to this class is 0%. 

• Class 1: one pair of lobes, with or without a torus-like structure. The fraction 
belonging to this class is 5(±4)%. 
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• Class 2: two pairs of lobes, with similar or different lengths. The fraction belonging 
to this class is 46 (±5)%. 

• Class 3: three pairs of lobes. The fraction belonging to this class is 49(±8)%. 

We have adopted some "usual customs" of morphological classification of PNe: if two pairs 
of lobes look nearly perpendicular and one of them is significantly shorter, then the short 
pair is usually identified as the torus; shorter lobes overlapped with longer ones are called 
inner lobes rather than extra pairs of lobes. The classification of each image has been 
done by four individuals independently, and the overall results are listed in Table HJ The 
percentages above are from the overall results, and errors correspond to variations among 
individuals, which are within expectation and actually reflect the real situation of previous 
classification works. 

This exercise shows that by assuming the number of pairs is at most three, only half of 
the PNe with three pairs of lobes will be correctly identified. The other half will have the 
number of lobes underestimated, and 5% of them will even be misclassified as single-paired, 
or more commonly called, bipolar. In other words, the true number of those with three 
pairs of lobes should be ~2 times the number of visually identified ones. Denoting n(a, b) as 
the number of PNe with a-pairs that are intrinsic but only with 6-pairs that are observed, 
the above statement can be expressed as 



where N is the maximum number of pairs in all PNe. This is a recurrence relation of x, to 
be solved with observational data and simulations together. Although observational images 



3 



n(3, m) ~ 2n(3, 3). 



m=0 



In general, the true number of those with x-pairs intrinsic should be 
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are available in many PN catalogues, the projection effect makes it hard to tell whether an 
object with apparently x-pairs of lobes should be put in the first or third term on the right 
hand side of the above equation. Therefore at the moment it would be unwise to form a 
conclusion on the real fraction of multipolar PNe. On the other hand, it is expected that 
the fraction of apparent multipolar PNe 

X I X 

^^n(x,m) / n(x, m) 

m=2 I m=0 

will increase with the actual number of pairs x, and thus the number of nebulae with over 
three lobes is more closely estimated than the three-paired ones. 



4.2. The Effect of Sensitivity 



In real observations, the morphological classification also suffers from sensitivity issues. 
In images of limited sensitivity or dynamic range, fainter features will be missed and their 
multipolar nature not recognized. Figure [2] illustrates a comparison between the images 
under high and low sensitivity conditions. The second row images are the same as the top 
row except that the pixels with brightness less than 1/3 of the peak brightness are not 
displayed. If one were to classify the apparent morphology based on the second row images, 
then one would arrive at very different morphological classes for the objects. 

Moreover, the "broken" segments of faint lobes create an illusion of the presence 
of minor structures. The bottom left panel in Figure |2] looks like a bipolar nebula with 
point-symm etric ansae, filaments o r knots around t he lobes, similar t o the Cat's Eye Nebula 



(NGC 6543, 



Balick fc Hajian 



2004J) and NGC 3242 ( IRuiz et al 



201ll ). In fact, these features 



are the overlapping areas of two or more lobes: overlapping increases the column densities 
of walls so that they can still be seen while the other parts of the lobes are filtered out. 
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4.3. Comparison with Observations 

In this section, we explore whether the simulated images are of any relevance to real 
observed PNe. We match the 20 selected real PNe images with simulated images of our 
simple model. While we do not claim that the sample is representative of all multipolar 
PNe, we hope to demonstrate that our model can be the first-order approximations to the 
true morphologies of the selected objects. For simplicity, we change only the six angles as 
described in Section 14.11 and keep all of the other parameters fixed. The effect of sensitivity 
is not included. Comparisons between the observed and simulated images are presented in 
Figure [3j The parameters used in the simulated images are listed in Table |2j The labels of 
the objects in Figure [3] are given in the first column of Table EJ 



4.4. Individual Objects 



Details of each object in Figure |3] are given below. Here, we try to relate the observed 
morphological features to the 3D model. 



(a) IRAS 05028+1038 



excitation PN ( IHarman et al 



and mass M C g = 0.79M© (McCarthy et al. 



320, PN G190.3-17.7) was classified as a Type II low- 



2004T) with central star effective temperature T e g = 85000K 



19971 ). It was describe d as a "poly-polar" PN 



200J). 



with surrounding knots distributed in a point-symmetric pattern (IHarman et al. 
The HST image clearly shows two knots both to the north and the south, aligned not 
exactly in the direction of the central parts. In our model we align the two pairs of lobes in 
the direction of the knots and align the remaining one to the central lobes in the NW— SE 
direction. Knots are regarded as higher order properties that can be represented by 
en hancing densities at t he tips, or as results due to sensitivity (see Section l4~2l) . According 



to 



Harman et al. 



(120041 ). there is one more knot further away to the west. This may be 
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related to our third pair. 



(b) IRAS 07172-2138 



a point-symmetric shape (ISahai et al. 



(Ml-12, PN G2 35.3-03.9) has closed collimated pairs of lobes in 



201lh . The central region is due to two overlapping 



pairs and the extended features in the NW-SE direction are explained as the third pair. 

(c) IRAS 10197-5750 (Hen 3-404, OH 284.2-0.8, Roberts 22) was described as a bipolar 
D UPLEX (DUst- Pro minent Longitudinally-EXtended) nebula with an A2 I central star 



by 



been dete cted ( 



Ueta et al.l ( 20071) . Fluorescence- induced H? emission, characteristic of young PNe, has 



Garcia- Hernandez et al. 



emission (Allen et al. 



20021 ) . It is also associated with intense OH maser 



19801 ). The equatorial dark lane in the optical image, together wit h 



the near- and far-infr ared (IR) image, implies the presence of a dusty ring ( 1 Allen et al.l 



1980 



Cox et al. 



20111 ). which is not included in our model. 



(d) IRAS 10214-6017 (Hen 2-47, PN G285.6-02.7) was classifi ed as a young, low- 



excitation PN with cool IR continua and [Ne n] 12.8 jum emission f lVolk fc Cohen 



1990). 



It has six obvious closed lobes which appear to be axis-symmetric about the equa torial 



direct ion. The components in the east side have a larger size than the oth ers, an d 



Sahai 



20001) is 



(120001 ) further divided them into two more lobes. The central bright ring (jSahail 
reproduced in the modeled image as the overlapping regions of the three pairs, assuming a 
simple addition of the column density. 

(e) IRAS 15015- 5459 (Hen 2-115 PN G321. 3+02.8) has a central star of T eff = 27400K 



and M C g = O.68M (IZhang fc Kwok 



19931 ) surrounded by warm dust. In addition to a pair 



of elon gated lobes, there are point-symmetric bulges nearer to the center. 



Sahai fc Trauger 



( 119981 ) identified a small rhomboidal structure around the central star, which may also be 



an illusion due to overlapping of the multiple pairs in our model. 



(f) IRAS 16409-1851 (Hen 2-180, PN G000. 1+17.2) shows two lobes slightly bent 
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to the west. Such "bending" structures are very difficult to understand by traditional 
physical ejection models, but can arise naturally in this simulation: the bent lobes seen 
at the west side belong to two pairs of symmetric lob es but their count erparts in the east 



20111 ). which is again a 



are significantly fainter. It has a narrow bright waist ( jSahai et al. 
possible consequence of superposition of multiple lobes. 

(g) IRAS 16585-2145 (IC 4634, PN G000.3+12.2) has a curious shape consisting of an 



S-shap e main nebula with extended S-shape bow-shock structures studied by 



Guerrero et al. 



( 120081 ). In the central region, an inner shell is hidd en roughly aligned wi th the main axis 



and kno t s are found along the bow-shock features ([Guerrero et al 



20081 ). According to 



Schwarzl ( 119931 ). the central part exhibits both red and blue shifted features on the same side 
along the main axis, suggesting the presence of a fast precessing collimated outflow with 
time-dependent velocity for shaping the nebula. The S-shape point-symmetric structure of 
the nebula can be explained by two overlapping bipolar lobes and the observed S-shape can 
be reproduced by non-uniform density distribution along the azimuthal angle of the lobes. 

(h) IRAS 17028-1004 (Butterfly Nebula, M 2-9, PN G010.8+18.0 ) has highly 



collimated open-ended lobes with skirt-like structures (jSahai et al. 



20111 ). Bright k nots, o r 



ansae , are found along the long axis, possibly formed at the head of an early jet (ISoker 



1990l ). The 



studied by 



proper motions of t he dusty blobs much further away from the main nebula are 



Corradi et al. 



( 120111 ). Their work postulated a rotating ionizing beam to explain 
the ionized gas emission phenomenon in this nebula, and supported the idea that the 
appearance of multiple lobes was the result of an excitation gradient. However, considering 
that other nebulae with similar apparent morphology may not have been studied in as much 
detail as M 2-9 to reveal the physical properties, we present the model here to suggest an 
alternate possible choice. The modeled image suggests that the fork-like ends of the main 
lobe can be geometrically reproduced by two pairs of lobes, and the ansae (knots) are at 



13 



the interacting points of the three pairs. There are fine structures such as arcs not included 
here. In addition, the lobes are more likely to be open-ended (see Section l5TTj) . 



(i) IRAS 17156-3135 (PN G354.5+03.3) In low dynamic range images, the nebula 
has the typical shape of a bipolar nebula. The HST image reveals another pair of bipolar 
lobes with an axis aligned at an angle with respect to the primary bipolar lobes. We put 
the third pair of bipolar lobes to be pole-on in order to produce the circular feature at the 
center. Given three pairs of lobes, it would not be uncommon that one pair would lie close 
to perpendicular to the plane of the sky. 



(j) I RAS 17296-3611 



lobe pair ( jSahai et al. 



PN G351.9-01.9) has an inner bubble inside the main bipolar 



201lh . Considering the brightness difference between the central part 



and the lobes, we put two pairs as nearly pole-on. 

(k) IRAS 17389-2409 (Hen 2-267, M 2-14, PN G003.6+03.1) is a point-symmetric PN 
with bulges. The apparent shape is similar to IRAS 15015-5459, but with a bright waist 
that is not perpendicular to the largest pair of lobes. 

(1) IRAS 17410-3405 (Hen 2-271, M 3-14, PN G355.4-02.4) is an asymmetric PN with 
apparently more bubbles in the S— E side than the N— W side. The Ha brightness peaks at 
the central region and drops slowly outward. The NW bubbles are fainter and not as clearly 
seen as the SE ones. This asymmetry can be due to non-uniform brightness distributions 
on the two sides, or can be the result of circumstellar dust extinction. Besides the three 
pairs, there is another small bubble in the equatorial direction. 

(m) IRAS 17496-2221 (Hen 2-299, M 1-31, PN G006.4+02.0) has a peculiar hexagonal 
shape that can be reproduced by three pairs of bipolar lobes, with a pair near the plane of 
the sky representing the long axis of the nebula. The bright waist represents the overlapping 
central parts of the three pairs of lobes. 
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(n) IRAS 17549-3347 (Hen 2-313, PN G357.1-04.7) The central part of this nebula has 
a circular shape which is considered as a pair of nearly pole-on lobes. The more extended 
structure has a pointed S-shape. 



(o) IRAS 17567-3849 (Hen 2-320, PN G352.9-07.5) has a bri ghter, shorter 



Dennis et al. 



pair of 



(|2009|). 



lobes and two narrower, longer pairs described as "shoulders" by 
These authors suggest that this structure can be produced by bow shocks in their 
magnetohydrodynamic nested-wind simulations. More extended collimated "nose cones" 
have also been found. Such double "inner" and "outer" bipolar lobes can be naturally 
explained by a multipolar structure as seen in the simulated image. 



(p) IRAS 18022-2822 (Hen 2-339, M 1-37, PN G002.6-03.4) has three pairs of lobes 
distributed in a point-symmetric manner. An arc-like feature to the north is believed to be 



part of a halo (ISahai et al. 



20111 ). The central part seems to be hollow and a bright star 



can be seen. This classic multipolar structure can be reproduced by three pairs of bipolar 
lobes almost on the plane of the sky. 

Dits faint extensions around 



Sahai et al. 



(120111 ) probably 



(q) IRAS 18039-2913 (Hen 2-346, PN G002. 1-04.2) exhi 
a bright irregular center. The close-ended barrel recorded by 
refers to the center. Three bullet-like features can be seen in the SW part and there are 
faint lobes in the opposite direction. The three lobes in our model image point to these 
features but fail to reproduce the central barrel. An extra component may be needed. 

(r) IRAS 18430-1430 (M 1-61, PN G019.4-05.3) has one pair of narrow lobes stretching 
out from the point-symmetric central part which is possibly another two pairs of lobes. 

(s) IRAS 19431+2112 (Hen 2-447, PN G057.9-01.5) has a bright waist and point- 
symmetric lobes. The sh ape of two lobes are m ore obvious in the NE direction. Knots and 



filaments were identified (jCuisinier et al 



20051 ). 
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(t) IRA S 20090+3715 (He n 2-456, NGC 6881, PN G074.5+02.1) has a complicated 



morphology. 



Kwok fc Sul (120051 ) found four pairs of rings in the multipolar lobes aligned 



with one pair of lobes , and three rings at the waist aligned to another lobe pair. 



Guerrero fc Manchadd (119981 ) identified a loop to the SE of the main nebula, almost 
perpendicular to the major axis. The H2 image shows a much more extended size than the 
optical one, and the lobe to the NW is ev en larger than its counterpart, turning into an 



irregular shape ( IRamos-Larios et al. 



20081 ). O n the other hand, the 6 cm radio morphology 



is more elongated in the equatorial direction ( lAaquist fe Kwok 



19901 ). and only images the 



region with the highest emission measure. Our model can simulate the multipolar lobes but 
not the central bright region. 



5. Discussions 

PNe have many different morphological features which are often attributed to different 
physical regions created by separate physical processes. Our models suggest that many of 
the observed features can be explained by a single, unified model. Of the bright torus seen in 
many PNe, some are due to real volume density enhancements, but some probably represent 
the overlapping region of different lobes. Simulations similar to the one presented here can 
sort out the "false torus" percentage in the whole PN population. The point-symmetric 
S-shape morphology can be explained by nearly aligned pairs of lobes. Ansae may be the 
tips of unseen lobes. The double inner— outer bipolar lobes can arise from overlapping 
multipolar lobes. Due to the long path length, a pair of near pole-on bipo lar lobes can 



2008|). The observed 



easily be mistaken for a torus, as in the case of NGC 6720 ( iKwok et al.l 
images of PNe can be successfully modeled with a very simple multipolar model. The 
model lobes have no density variation (either radial or azimuthal) and all three lobes have 
identical shapes and sizes. There is no inclusion of the effects from dust extinction, which 
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can make the front lobes brighter and the back lobes fainter. 



5.1. The Shape of the Lobes 

The shapes of the lobes reflect the dynamical history of the outflows and the 
interactions among outflows. Assuming a collimated fast outflow from the central star, 
the shapes of the lobes are defined by the interaction between this fast outflow and the 
ambient circumstellar materials, e.g. the remnant of the AGB wind. At a later stage of 
development, the fast outflow may break out of the circumstellar envelope and change into 
an open-ended shape. 

The lobes are set to be close-ended in our model. The difference between close- and 
open-ended lobes is easily identified even in a projected image, unless the lobe is tilted 
almost pole-on. This can be seen by comparing Figure [T] and Figure HI Moreover, the 
identification of multiple lobes does not depend on this factor. Figure [5] shows another 
example of the closed-lobed model with a different lobe shape for comparison, described 
by the equation r = cos 9 in spherical coordinates. Compared to Figure [U Figure [5] shows 
greater deviations from the real PN images. In fact, lobes with larger opening angles have 
a higher chance of interacting with each other and producing more complicated shapes. 



5.2. The Number of Pairs 



In general, "multipolar" means having more than one pair of lobes, and is not confined 
to objects with three pairs. For the objects chosen, at least three pairs are obviously 
seen, and the number three is also c ommonly found i n the literature (e. g., NGC 7027 by 



Nakashima et al. 



2010 



NGC 6644 by 



Hsia et al. 



2010 



and NGC 7026 by 



Clark et al. 



It is possible that there are more than three pairs (e.g. IRAS 19024+0044 by 



2012 



Sahai et al. 
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2005 



and NGC 5189 by 



Sabin et al.ll2012l ). but adding more pairs means adding more 
parameters; at this stage we hope to keep the number of parameters to a minimum. The 
fourth pair is usually not as obvious as the other three. For IRAS 19024+0044, the fourth 
pair almost overlaps with another one, which could be due to projection effects, or the two 
components may have been in the same lobe produced by perturbations. At least three 
is enough to describe such objects. The less obvious lobes are therefore treated as higher 
ordered structures. 



5.3. Lobes: Same Length or Not? 

By looking at a 2D image, it is hard to tell whether the lobes have the same length or 
not without knowing the inclination angles. One of the possible ways to predict the length 
ratios is to compare the brightness of various lobes. If the brightness is simply proportional 
to the square of the column density, and if the volume density is constant everywhere, 
then the lengths along the LOS can be calculated. As shown in the modeled images, lobes 
that are more pole-on look brighter. However, inf rared images con firmed the existence of 



denser tori in som e pPNe, e.g., IRAS 16594-4656 (jVolk et al. 



flVolk et al. 



20061 ) and IRAS 17441-2411 



20071 ). which may be optically thick in the visible light. Even if the central 
bright ionized region is optically thin in PNe, it is still hard to tell whether it is a dense, 
short component or a less dense but longer one. 

Another way to determine the actual length ratios is from kinematic information. 
Lobes with the same lengths are likely to be produced simultaneously and with the same 
outflow velocities. Assuming constant expansion velocities, the ratio of the apparent size 
to the velocity component along the LOS should be able to indicate the projected angles 
of each lobe. Take i to be the angle between the tilted lobe and the LOS. The apparent 
angular size ?" a pp is related to the real size ?" rea j by a factor of sinz, while the LOS velocity 
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and real velocity are linked by cosi: 

rapp = r real sinz, (2) 

v LOS = w real cosz - ( 3 ) 

To find r rea Y and f rea ] one must know i. For one single pair of lobes, it is hard to deduce 
i, but with multiple lobes, the ratios between two pairs can eliminate the real quantities, 
based on the assumption of the and the same v Teg ^ for all lobes. On the other 

hand, for lobes of different lengths, deducing the true length ratios from 2D images is more 
uncertain. 

A detailed morpho-dynamical study is the key to a reliable determination of the real 
shape of a PN. The kinematic information is the only way to "see" the third dimension 
in the 2D sky plane. In addition to the estimation of the length of the lobes, visually 
overlapping lobes can also be distinguished from their radial velocities because they are 
tilted with different angles to the LOS. 



5.4. Formation and Evolution of Multipolar Lobes 



Although the formation mechanism of multipolar lobes is not confirmed, the theory 
of a precessing jet is more frequently adapted than coeval mass outflows. As in the case 
of bipolar PNe, a binary system is believed to be a driving source. The jet changes 
direction due to the a dditional velocity component of the progenitor's orbital motion 

and references therein) and describes a s piral trajectory, which can 



(IVelazquez et al. 



2012 



be traced by motions of water maser spots in an earlier phase ( jlmai et al. 



20021 ). It has 



been suggested that the number of lobes is proportional to the ratio of the jet precession 
period t o the orbital period, w hich depends on the masses of the two objects in the binary 



system (IVelazquez et al. 



20121 ). In this hypothesis, the angles between adjacent lobes should 
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be equal, and all the lobes should be confined to a cone swept out by the jet. Compared 
with our random direction simulation, this hypothesis puts more constraints on the lobe 
orientations, and therefore projection effects are expected to be more significant: the lobes 
have a smaller angular distance to each other and thus a higher chance of overlap in the 
projected images. Furthermore, if the precessing jet mechanism is true, then the S-shape 
structure found in some PNe is just a variation of the multipolar scenario, where the "lobes" 
are so close to each other that a continuous S-shape is seen instead of discrete lobes. 

In each lobe, expansion components in the direction perpendicular to the lobe axis push 
out the lobe walls. Eventually, the lobes interact with each other, combine into a larger 
lobe, or become less obvious because of the fusion of the walls. As a result, the phenomenon 
of multipolar lobes is more commonly observed in pPNe or young PNe. If the existence 
of a companion is the only way to produce multipolar outflows, then the total fraction 
of multipolar PNe can be more easily estimated as the product of the fraction of binary 
systems in all PN progenitors and the fraction of multipolar outflows onset in these binary 
systems, scaled by the ratio of survival duration of multiple lobes to the lifetime of a PN. 
Although the fraction of binary systems directly confirmed in individual PN progenitors 
is lower than the percentage of multipolar PNe (12%— 20%, see the introduction section), 



the overall estimated fraction o f binary progenitors is comparable: t he close binary 



might be as high as 12%— 21% ( IMiszalski et al 



20091 ) or 22% ± 9% flFrew fc Parker 



r action 



201fl) 



6. Summary 

From a simple three-bipolar-lobes model, we have shown that many of the observed 
morphological features of PNe can be successfully reproduced. Morphological features such 
as "point-symmetric" , "S-shape" , "bending" , "double bipolar lobes" , "knots" , and "ansae" 
can arise naturally from this model. These simulations confirm the importance of the 
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effects of orientation and sensitivity, as discussed in iKwokl ( 120101 ) . We have also shown that 
the statistical analysis of morphological classes based on apparent shapes is not a reliable 
indicator of the true structures of the nebulae. 

The reconstruction of the true 3D structure of PNe represents the first step in the 
identification of the physical processes responsible for the shaping of PNe. Without an 
accurate picture of the 3D structure, all attempts at physical models are futile. 

Further expansion of the model should take into account the density profiles, ionization 
and radiative transfer treatments, dust extinction, and interaction with the circumstellar 
and interstellar media. The fourth dimension — temporal evolution is just as important. 
The 20 PNe presented are barely the tip of the iceberg among the known multipolar PNe, 
and our simulations suggest that the number of multipolar PNe will keep on rising as the 
sensitivity and dynamic range of images improve with future observations. 



We thank Nico Koning and Wolfgang Steffen for discussions on the use of the SHAPE 
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v-— / s + \ x * k * 

Fig. 1. — 100 projected images of the three-lobed model with randomly generated inclination 
and position angles for each pair. The scales are the same in all panels. 
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Fig. 2. — Perception of morphology is affected by viewing angles and sensitivity. Upper 
row: with the six angles in the model fixed, the viewing angle is changed from each image 
to the next by 15° of % and 15° of PA together. Lower row: each image is modified from 
the one above so that the faintest pixels below one-third of the peak brightness are cut off. 
Brightness levels are shown in linear scale. 
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Fig. 3. — Comparison of the 20 observed images (those with an index letter), each with its 
corresponding modeled image (the one right below each observed image). Brightness levels 
are in log scale. North is pointing up and east to the left. Refer to Table [2] for the object 
data. The angular sizes and brightness levels of the objects are not the same. 
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Fig. 5. — Same as Figure [T]but with a different lobe shape. See the text for the descriptions. 
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Table 1. Randomly Generated Angles of Each Image in Figured], from Left to Right, from 

Top to Bottom. 
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Table 2. Information and parameters of the 20 PNe. 
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-84 


24 


7 


28.2 


29.6 


46.2 


1 


17410-3405 


U6MG3101B 


25 


-40 


32 


-23 


25 


12 


10.5 


17.7 


21.6 


m 


17496-2221 


x TrTTTT/^Ann x~* 

U5HH6902B 


40 


3 


23 


-35 


27 


33 


20.6 


25.7 


27.3 


n 


17549-3347 


U6MG3601B 








31 


18 


22 


-14 


16.6 


21.6 


31.2 





17567-3849 


U5HH1302B 


41 


82 


35 


71 


10 


79 


8.5 


25.1 


30.3 


P 


18022-2822 


U35T2105B 


37 


-12 


39 


-50 


33 


-77 


17.5 


23.3 


36.1 


q 


18039-2913 


U5HH4103B 


41 


-3 


45 


-26 


58 


-42 


16.0 


18.4 


34.0 


r 


18430-1430 


U59B0301B 


21 


-62 


52 


-19 


20 





20.5 


33.2 


38.4 


s 


19431+2112 


U59B0704B 


51 


44 


50 


67 


27 


65 


17.8 


23.5 


27.7 


t 


20090+3715 


U39H3601B 


50 


-22 


40 


-55 


80 


-33 


25.1 


31.6 


43.9 



-36 - 



Note. - - The order of pairs 1, 2 and 3 is arbitrary. 9 (minimum, median or maximum) 
refer to the separation angles between pairs. 



